Stress-Energy Tensor Induced by Bulk Dirac Spinor in Randall-Sundrum Model 
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Motivated by the possible extension into a supersymmetric Randall-Sundrum (RS) model, we 
investigate the properties of the vacuum expectation value (VEV) of the stress-energy tensor for 
a quantized bulk Dirac spinor field in the RS geometry and compare it with that for a real scalar 
field. This is carried out via the Green function method based on first principles without invoking 
the degeneracy factor, whose validity in a warp geometry is a priori unassured. In addition, we 
investigate the local behavior of the Casimir energy near the two branes. One salient feature we 
found is that the surface divergences near the two branes have opposite signs. We argue that this is 
a generic feature of the fermionic Casimir energy density due to its parity transformation in the fifth 
dimension. Furthermore, we investigate the self-consistency of the RS metric under the quantum 
correction due to the stress-energy tensor. It is shown that the VEV of the stress-energy tensor and 
the classical one become comparable near the visible brane if k ~ M ~ Mpi (the requirement of no 
hierarchy problem) , where k is the curvature of the RS warped geometry and M the 5-dimensional 
Planck mass. In that case the self-consistency of RS model that includes bulk fields is in doubt. If, 
however, k < M, then an approximate self-consistency of the RS-type metric may still be satisfied. 

PACS numbers: 04.50.-h, 04.62. +v, ll.10.Kk, ll.25.-w, 11.30.Pb 



I. INTRODUCTION 

The hierarchy between the Planck scale, Mpi ~ 10 19 
GeV, and the standard model (SM) scale, Msm ~ 1 
TeV, has been a long-standing problem in high energy 
physics. In the past decade, there have been two popular 
solutions to the hierarchy problems: the Arkani-Hamed- 
Dimopoulos-Dvali (ADD) model 1 and the Randall- 
Sundrum (RS) model [2J. Both models invoke extra di- 
mensions and the brane- world scenario. The weakness 
of gravity is associated with the largeness of the extra 
dimension in the case of ADD, while in RS it is due to 
the exponential warpage of the extra dimension. In this 
paper we shall only focus our attention on the latter. In 
the RS model, two flat, parallel 3-branes are located at 
two fixed points on a S 1 /Z 2 orbifold. The metric of RS 
reads: 

ds 2 = e'^^rj^dx^dx" + dy 2 . (1) 

where R is the radius of the orbifold and y is the extra di- 
mension coordinate ranging from — nR to -kR and <j(y) — 
k\y\. We adopt the convention 77^, — diag(— 1, 1, 1, 1). 
The brane located at y = is called the hidden brane 
and that at y — irR is the visible brane, on which the 
familiar SM fields reside. Under this construction, the 
hierarchy problem is naturally solved via the warp factor 
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e -kixii a i on g the extra dimension, generating a large hier- 
archy without requiring large extra dimensions. Specif- 
ically, with the choice of kR ~ 12, the TeV scale at the 
visible brane can be descended from the Planck scale at 
the hidden brane. 

With the introduction of a finite separation between 
the two branes in the extra dimension, any field in the 
bulk should induce a Casimir energy. This provides a 
possible resolution to the smallness problem of the cos- 
mological constant when connected with the observed 
dark energy [3H5] ■ The Casimir effect in RS model gener- 
ated by scalar fields has been studied in details by many 
authors [6Tfl3]. In [6], in particular, a detailed analysis 
of the local Casimir energy, or to be more specific, the 
vacuum expectation value (VEV) of stress-energy ten- 
sor, is given, in order to test the self-consistency of the 
model. It should be noted, as was pointed out by the au- 
thors of [BJ, however, that some conclusions made therein 
might not hold in a supersymmetric RS model, where 
the fermionic contribution is expected to exactly cancel 
its superpartner counterpart yet its local behavior has 
never been explicitly investigated. It is thus desirable to 
investigate the effect of the Casimir energy induced by 
the fermion field in the bulk. 

With regard to the fermionic Casimir effect in RS ge- 
ometry, some authors obtained the result by multiplying 
a degeneracy factor to that of the scalar field [HI [16] . 
In such an approach, calculations are performed via the 
mode summation method. We note that while such a 
transcription from the scalar field Casimir energy to that 
of the spinor field in a flat spacetime geometry is clear, 
that for the curved spacetime is not as transparent. This 
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is mainly because the reduction from the Dirac equation 
to the Klein-Gordon equation in flat spacetime is un- 
ambiguous, whereas the similar procedure in the curved 
spacetime would induce a curvature coupling term |15j : 



{-V 2 +m% + ^R)^ = Q, 



(2) 



where ip is the fermion field, R is the Ricci scalar and 
to,/, is the mass of the fermion. On the other hand, for 
a scalar field <j) in the curved spacetime, the equation of 
motion reads: 



£R)4> = 0. 



(3) 



The ambiguity arises because the coupling parameter £ 
between the scalar field and the curvature is generally 
not restricted to but as a free input to the theory. 
This motivates us to calculate the Casimir energy via 
the Green function method where the stress-energy ten- 
sor can be derived explicitly without resorting to the de- 
generacy factor. 

We note that the vacuum energy arising from the 
fermion field in RS model has been fully discussed by 
Flachi et al. [17] [18] by computing the lowest order quan- 
tum corrections to the effective action. The finite tem- 
perature effect of the Casimir energy arising from both 
scalar field and spinor is discussed in [15] . These ap- 
proaches can only provide the global properties of the 
vacuum energy and important local behavior might be 
overlooked. With the interest to fully understand the 
nature of the Casimir energy in the RS brane world, we 
find it desirable to directly compute the VEV for the 
stress-energy tensor, so as to compare it with that in- 
duced by the scalar field explicitly in the context of a 
SUSY system. 



II. BASIC SETUP 

Let us first distinguish two sets of labels in our metric: 

ds 2 = g M Ndx M dx N = e-^^v^dx^dx" + dy 2 , (4) 

where the capital Latin labels denote the (4+1)- 
dimensional quantities: M — 0, 1,2,3,5, so that x° = y, 
whereas fi — 0,1,. .,3. Then we write down the action 
for a Dirac spinor in the bulk with mass mq, — c<j' 



d 4 x 



(5) 



where g — det((?MA r )- The gamma matrices, 7 are 
defined in curved spacetime as j M = e Q , M 7 a , where 
e a = diag(e CT , e CT , e CT , e a , 1) is the inverse vierbein and 
7™ = (7°, 17 s ) are the gamma matrices in flat spacetime, 
and Dm = 9m + Tm is the covariant derivative in curved 
spacetime. With the metric defined in Q, we have 



0. 



(6) 



It should be noted that the equation of motion derived 
from ([5]) with to* = is automatically conformally 
covariant under proper choice of boundary condition, 
without the need for an additional conformal coupling 
term [2"4]. 

Given the Z2 transformation for Dirac fermion: 
y) = ±J 5 ^(y), it can be seen that ^f^> is odd un- 
der Z2 transformation. This implies that to* must also 
be odd under Zi transformation in order to preserve the 
Zi symmetry of the Dirac equation. Therefore, to* can 
be parametrized as [30] 



771,5 = ca 1 = cfce(y). 



(7) 



where e(y) is defined as being 1(-1) for positive(negative) 

y- 

The stress-energy tensor related to the above action is 
given by 



Tmn = i^Af-tV)* - 9m nL, 



(8) 



where the second term that involves the Lagrangian does 
not contribute to VEV by virtue of the equation of mo- 
tion, and will be neglected in the following calculation. 



III. COMPUTATION OF THE GREEN 
FUNCTION 

In order to obtain the VEV of stress-energy tensor, 
we first calculate the Green function G(x M , x' M ) for the 
field, then express {Tmn) in terms of the Green func- 
tion. The Green function by definition satisfies the field 
equation 



{ij M D M -m*)G(x M ,x ,M ) 



M iM\ 



1 



S(x M -x ,M ). (9) 



To eliminate the dependence on the coordinates of x^ , we 
perform a 4-dimensional Fourier transform on the Green 
function 



G(x M ,x' M ) 



d 3 p 
(2tt) 3 



duj 
2^ 



e ip-(£-*') e -M*-t')G p (y )2 /). 

(10) 



Then ^ becomes 
[-7% - 7 5 3 5 - 771* ]G p (y, y') = e 2 °6(y - t/)/ 4x4 . (11) 



where G p (y,y') = e - 2 °G p (y,y'). 

We further write the Green function as 



G P (y,y') 



h h 

91 92 



(12) 



where fx, /2, 51,32 ar e 4x4 matrices. Put this back into 
(111, we arrive at the following 4 equations for the 4 
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elements of the Green function: 

-dsfi + ckh + e°a a p a9l = -e 2a 5(y - y')I 2x2 , 
e a v a p a fi + d 5 gi + ckg 1 = 0, 
-#5/2 + cfc/ 2 + e a a a p a g 2 = 0, 
e°v a p a fz + d 5 g 2 + ckg 2 = -e 2a 5(y - y')I 2x2 . 



(13a) 
(13b) 
(13c) 
(13d) 



If we neglect the Dirac delta functions on the right hand 
sides of (13a|~ (13d I, the solutions are [20] 



My) 
9i(y) 



[J c _i (mz) + b l H 1 c _ 1 _(mz)], (14) 
[J c+h {mz)+d l H l k {mz)]. (15) 



where m = \J—p 2 is the Kaluza-Klein mass, z = e CT /fc, 
i = 1,2. Ni and bi are the coefficients to be determined 
by normalization and by our choice of the boundary con- 
dition (BC). Generally, there are two distinct classes of 
BCs for Dirac fermion, the untwisted and the twisted 
BC[H]EI]- These BCs are derfined as follows. 



Untwisted: 

|o=0, 

(d 5 +ck)^ R | = 0, 
Twisted: 

*L | =0, 

{d 5 +ck)^ R | =0, 



*L U= 0. 

(d 5 + ck)t> R U= 0. (16a) 

{d 5 - cfc)* L \ nR = 0. 

U= 0, (16b) 



where ^> L R = [(1 =p 7 5 )/2] v I / . We'll deal with both cases. 



A. The Untwisted BC 

For later convenience, we introduce the following four 
functions to represent different combinations of the spe- 
cial functions that appeared in ( 14 ) or (fT5l 



viv) 

m = 



e 2 [ J c _ 1 (mz) + bH\_ 1 (mz)] , 

2 c 2 

[ J c _ 1 (mz) + bH\_ 1 (mz)] , 
-e% [J c+ i(mz) + bHl + i(mz)], 
-e f [^iW + bH l c+l (mz)], 



where 



b = 



(17) 
(18) 
(19) 
(20) 

(21) 
(22) 



In terms of these new functions, the untwisted BCs be- 
come 



V lo= 0, 
(d 5 +ck)X | =0, 



V \itR= 0. 

(*l#U=0. (23) 



We first deal with /1 and g\. Under the above definitions, 
fi and (71 can be expressed as 



h{y,y') = 





?) 






1(2 









(24) 



gi(y,y') = < 



v(y), if y > y' 



I JjA(y), if 2/ < 



/3 I 



(25) 



A(y), if y > y' 



All the parameters in the above matrices are determined 
by (13a) and (13b I, from which we get 



gi \y=y>+ ~gi \y=y>-= 0, (26) 
fl \y=V+ -fl \y=y'-=e 2 ^y'h 2x2 , (27) 

- dufi + ckh + e a a a Pa9l = 0. (28) 

These together constitute 16 linear equations for 
16 unknowns. Here we solved it by the help of 
Mathematica 7.0: 



h(y,y') 
9i(y>y') 

where 



-I 2x2 e 2 ^^f V (y), iiy<y> 
-I 2x2 e 2 °^^fri(y), if y > y> 

-a<*p a e 2 <y'^\(y), if y < y' 
-a<*p a e^v')^\(y), if y > y> 

S = Xl] - Xlj \y> . 



(29) 
(30) 

(31) 



The remaining two elements, f 2 and g 2 , can be obtained 
in a similar way 



, , % f -^P a e 2a{v,) ^Q-r,{y), if y < y' , 

fo(v v ) = < (32) 

, / -l2, 2 e 2 ^^fx(y), i£ V <v> 

fl2( ^ ) = \-w^)^),if y >^ (33) 



With (29), (30), (32), (33), we finish our calculation for 



the Green function subjected to the untwisted BC. 



The Twisted BC 



In order to fit the twisted BC (16b), we define addi- 
tional two functions 

f)(y) = -e^[J c+ i{mz) + bHl +1 _(mz)], (34) 

X(y) = e%[J c _i(mz) +bHl_ k (mz)], (35) 
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where 



T I ' 



(36) 



so that 



V \irR= 0, (<9 5 - ck)X | 



7T.R- 



0. 



(37) 



If we rewrite ( 24 ) and ( 25 ) with the replacements rj 



A and A — > fj, the remaining calculations are completely 
identical to that of the untwisted case. The results are 



just the same as (29), (30), (32) and (33) with the above 
substitution. 



IV. STRESS-ENERGY TENSOR 

Now well-equipped with the exact form of the Green 
functions for both untwisted and twisted BCs, namely, 
(33), together with the definition (12), we are fi- 



(29 



y at a position to calculate the VEV of the stress- 



nal 

energy tensor. First, we make the identification 

iG(x M ,x' M )) = (V(x M )*(x' M )). 



(38) 



With this identification, we can replace any VEV involv- 
ing the quadratic terms of the field by the Green function. 
We begin with the VEV of T 00 in ([§, 



(Too) =ie 00 (^°D y) 

= -ie- a [id Tr( 1 °G)] x M =x ,M. 



(39) 



In the second equality we have neglected the spin connec- 
tion term To in the covariant derivative Dq, since it turns 
out to be a term independent of the mode, and thus can 
be omitted in the renormalization. Let us first deal with 
the untwisted BC. We perform a 4-dimensional Fourier 



transform on both sides of (39) and use (12) to obtain 



(*0o) = -ioje a Tr(g 1 + / 2 )| y =y' 



-2ie 



3a 



U! 2 AA + rjij 

m Xrj- Xrj 



(Untwisted BC), (40) 



where t o is the 4-dimensional Fourier transformation of 
Too- In a similar fashion, we obtain the other entries of 
the stress-energy and also those for the twisted BC as 
follows. 



-2z^[^/f^g±M, untwisted 



(Tyy) — 

-^1 (0/ t^(l^S) 1,=,', untwisted 
-^/(0/^5 5 (g±^)|^, twisted 



(41) 



where <5 M „ = diag(l, 1, 1, 1). This is our final expression 
for the stress-energy tensor. The integral cannot be ex- 
pressed analytically for a general c, and suitable regu- 
larization must be implemented. For the massless case, 
however, the results can be written in a concise way. 



V. MASSLESS DIRAC FERMION 

Now consider the simplest case: c = 0, in which all rj 
and A functions reduce to basic trigonometric functions, 
and (41) and (42) become 



1)], untwisted 

2i<W 3CT f(f?r J^^F tan[f (e*** - 1)], twisted 

(43) 



(Tyy) — 



2^/(0/f mcot[f( e 

d 3 p 



m ( a -ivkR 



1)], untwisted 
1)], twisted 

(44) 

The integrations are best done by performing a Wick 
rotation 



w — > ip~, m = \fuj* — p* — > ip 
Thus the stress-energy tensor is given by 

r^^/^^cothiK^-i)], 



(45) 



-2^/^^tanh[£(e 



\V( e KkR _ 



untwisted 
1)], twisted 

(46) 



{ 2e 5 ^ I ^ I pcoth[f( e ^- 
be^/^ptanhlKe^ 



1)], 
1)1, 



untwisted 
twisted 



(47) 



The regularization procedure for the integral is conven- 
tional, where one subtracts 1 from the coth or tanh func- 
tion. The results read (for a more detailed calculation, 
see, for example, [6] or [25] ) 



\^jj,v/ rcn 



-e 3<T r^2- 3 7r-f a- 5 r(f)C(5) , untwisted 



15 3<r 

16 c 

(^2/2/)ren 



?^2- 3 7T-5a- 5 r(|)C(5) , twisted 



ie 5ff 7r-§a- 5 r(f)C(5) 
-||e 5CT 7r-fa- 5 r(|)C(5) 



untwisted 
twisted 



(48) 



(49) 



where 



( „-irkR 



(e 



l)/k. Compare with the VEV of 



the stress-energy tensor for scalar fields given in [6], we 
conclude that 

Dirac, untwisted — real scalar 

(50) 

15 



AIN) Dirac. twisted 



-(Tmn) 



real scalar 



(51) 
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The factor 4, as indicated in [H| . accounts for the differ- 
ence of the degrees of freedom between the Dirac spinor 
and the real scalar fields. The sign difference originates 
from the distinct natures of fermions and bosons. As for 
the factor 15/16 between the untwisted and the twisted 
results, it results from the difference in regularization be- 
tween the coth and the sinh functions. 



VI. MASSIVE DIRAC FERMION 

To extract useful information from the VEV of stress- 
energy tensor for a general mass c, we need to appeal 
to the numerical method. In this section, we will only 
focus on the numerical integration of the 00 component 
of (Tmn) for the untwisted BC (the general features, 
including the power and the signs of surface divergences, 
are the same for the twisted BC) , namely, ( 40 ) . For later 
convenience, we denote 



F( P ,y,R) = 



XX 



r\r\ 



Xr] — A// 



(52) 



Recall that R is the radius of the RS geometry. So our 
goal is to regularize the integral 



(T 00 )=iC dpp 4 F(ip,y;R), 
Jo 



(53) 



0.000015 - 



0.00001 - 



5.x 10~' 



0.5 



1.5 

y 



2.5 



FIG. 1: The massless (c—0), untwisted BC case for (Tbo)rcn 
in units of k 5 . The branes are located at y = and y — n 
and kR = 1. Note that the numerical value (points) agrees 
perfectly with the exact solution. 



o.oooi P 



-0.00005 



-0.0001 - 



0.5 



with 



C 



r(§) 



4 r(3) 



=>3<r(l/) 



(54) 



FIG. 2: The massive (c = 10~ 3 ), untwisted BC case for 
(2oo) r en in units of k 5 . The energy density diverges at both 
branes, but with opposite signs. 



Note that the Wick rotation has been performed in the 
above expression. 

To regularize the integral, we subtract the same quan- 
tity, namely, 1, as the massless case from the integrand 
iF(ip,y; R). In other word, since iF{ip,y\ R) approches 
unity as p — > oo, we are in fact subtracting its ultraviolet 
limit [27]. 



C dpp 4 [iF(ip,y;R)-l] (55) 



(Too) i 



To facilitate the computation, we expand iF(ip 1 y; R) 
to the hrst order in c. Since c represents the fermion mass 
in units of Planck mass, it is expected to be very small 
for a physical particle, thus ensuring the validity of our 
approximation. In addition, to make the computer run 
more effectively, we set k = R = 1 as in [B], although the 
more "realistic" value would be kR ~ 12 with k ~ Mp\. 

The results are given in Fig. [I] and Fig. [2] It can be 
seen that the energy density diverges to infinity near both 
branes, but with opposite signs. We'll discuss this feature 
in the following subsections. 



A. Surface Divergence 

In Ref.[2S], the authors pointed out that the surface 
divergence is mainly contributed by the high wave num- 
ber (momentum) modes in the Fourier transformation of 
(Tmn}- Base d o n the same philosophy, we expand the 
integrand of ( 55 ) asymptotically in powers of p~ n and 
retain only the p~ l term, which dominates the integrand 
for large p, i.e., 



POO 

(T o) rc „ ^ C J dpp 4 [coth(— ^— )p - 1 



CI dpp 4 -^— 
Jo 2za P 



csch( 



a - 1 



)P 



x |z(a — 1) + acosh[2p(y — z)] — cvcosh[2p(z — 

(56) 

where a = e 7vkR , and recall that z = e a jk. The first term 
is just the result for massless fermion while the second 
term corresponds to the correction due to the fermion 
mass. Note that the mass correction term is proportional 
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to c by virtue of the small mass expansion to the first 
order. 

The expression remains complicated. However, since 
only the high wave number behavior concerns us, we iso- 
late the dominating term by taking the large p limit again 
and find 



00/SD 



C 



dp 



2cp 



,3 r 



-2p(z- 



-2p(f-z) 



(57) 



(The subscipt SD stands for the surface divergent term) . 
Since z = ranges from \ to a = ^j— , we find the 
near-brane behaviors to be 



ce 3ky 

(^bo)rea|y->0 ^ 4 : 



(Too) Ten | y— >7tB. ^ 



- ce 3ky 



(58) 
(59) 



We observe that while the energy density diverges to mi- 
nus infinity at the visible brane, that near the hidden 
brane approaches to plus infinity. This surface divergence 
agrees qualitatively with our numerical result shown in 
Fig. [2] It should be noted that the divergence only disap- 
pears in the limit of c = 0, the massless and conformally 
symmetric case, which is consistent with our previous 
result. 

At the first glance, the asymmetry aspect of the surface 
divergence seems curious. One might blame the asymme- 
try to the nature of the RS geometry. However, this is 
not the case; the fermionic Casimir energy density in the 
flat case k = also possesses such asymmetric surfact 
divergence. 



B. The k = Case 

In the k = flat case with untwisted BC, where we set 
the two branes to be located at y = and y = oi, the 
(17l~(20l should be replaced by 



n(y) = smK(y-a ), 
V(v) = sinK(y - eii), 



to* 



\{y) = cosK(y - o ) — sinK(y - a ), 

A 



X(y) = cosK(y - a x ) 



m* 

Ik 



smK(y -ai), 



(60) 
(61) 

(62) 
(63) 



where K = -J to 2 — to 2 , and m is the KK mass defined as 



before. Putting the above functions into (52) and (53), 
performing a Wick rotation upon to, we find a neat ex- 
pression for the k = case. 

(Too) fe =o = C(k = 0) f°° dpp A (TW+TM+T&>), (64) 



where 
T (o) 

T (i) 

T (2) 

with k 



coth[«(ai - a )y], 

sinh[«;(2y — ai 



a )} 



k sinh[K(ai — oo)] 
/to* \ 2 smh[n(y — ao)]sinh[K(y — oi) 
V K ) sinh[Ac(ai — oq)] 



(65) 
(66) 

(67) 



Vp 1 



? Tr , the Wick rotated variable of K. 



The superscript (i) stands for the ith order term (if not 
counting the k dependence) of the fermion mass to*. 
Taking the ultraviolet limit, we find that the latter two 
terms contribute to the surface divergence feature: 



00/fc=0,SD 



dp p 



4 m* 



e -2«(v-oo) g-2«(ai-») 



+ 



p 1 /TO* 

dp Y{— 



e -2K(y-a ) , e -2«(ai-j/) _ ^ 



(68) 



The first line comes from T*- 1 ', the first order mass term, 
and the second line comes from T< 2 ). The minus one m 
last parenthesis is independent of a\ and ao, thus could 
be eliminated in the rcgularization. Compare the lead- 
ing order term in (68) with (57), we note that the surface 



divergence features are basically the same, so the asym- 
metry aspect doesn't seem to originate from the nature 
of RS geometry, but a generic feature for a 5D fermionic 
Casimir energy density. In addition, we note an inter- 
esting fact in (68): the surface divergent term is anti- 



symmetric with respect to the two branes for the first 
order mass term, and symmetric for the second order 
mass term. 

In fact, the asymmetry aspect of the energy density 
originates from the parity transformation y — > —y. Un- 
der the parity, the fermion field mass changes sign, 



m* — > —to*. 



(69) 



As a result, all the terms proportional to the odd power 
of mass will change sign, while those proportional to the 
even power remain the same under the parity transfor- 
mation. Thus the terms in the energy density will be 
antisymmetric for odd power of to*, and symmetric for 
even power ones. This is in accordance with our result 
(68) and (57). In other words, the sign of the mass will 



determine a specific direction in the Bfth dimension and 
thus break the exchange symmetry of the two branes. It 
is also worth noting that such asymmetry aspect doesn't 
arise in the scalar field [H[57], since the mass of the scalar 
field maintains its sign under parity transformation. 

With regard to the non-integrable divergence of the 
stress-energy tensor, the problem was solved by Kennedy, 
Critchley, and Dowker for a scalar field in a static 
spacetime, and further elaborated by Romeo and Sahar- 
ian [30J [SI] , and put in a broader context by Fulling [32] . 
The resolution lies in the renormalization of the bare sur- 
face gravitational action, which induces a delta function 
and that cancels exactly with the surface divergence. 
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Since the surface divergence should be renormalized 
into the surface action terms, we do not expect the pow- 
ers of surface divergence for scalar field and Dirac spinor 
field to be the same, so that they could cancel each other 
in a supersymmetric theory. In fact, while the stress- 
energy tensor for a massive scalar field diverges cubically 
near the surface (see (2.39a) in [27]), that for a massive 
Dirac spinor field diverges quartically as in ( 58 ) and (p39| , 



and that for a minimally-coupled massless scalar field 
diverges as inverse fifth power of the distance (see [6]). 
The differences could be understood by simple power- 
counting analysis. 

Generally speaking, a surface divergence term in flat 
(4+l)-dimension takes the form 



(T , 



dpp A —u(py), 



(70) 



with 17 the quantity, the mass for example, that brings 
about the surface divergence, and 7 the mass dimension 
of it. Since the Casimir energy density has mass dimen- 
sion 5 in (4+l)-dimension, u(py) is dimensionless. Con- 
sider a minimally-coupled massless scalar field, the quan- 
tity signaling the surface divergence is the dimensionless 
coupling constant, which means that f2 is dimensionless 
and 7 = 0, so the integration yields a y~ 5 surface diver- 
gence. On the other hand, f2 = ck = for a massive 
Dirac spinor field, so 7 = 1 and the surface divergence 
is of inverse fourth power of the distance from the sur- 
face. As for the massive scalar case, since the mass of a 
scalar field generally appears in power of m^, Q = 
and 7 = 2 for a massive scalar field, hence the surface 
divergent term behaves as y~ 3 . 



VII. DISCUSSIONS AND CONCLUSIONS 

In flat spacetime, the one-loop correction to the ef- 
fective action, and therefore the Casimir energy, of the 
massless Dirac spinor field can be obtained simply by 
multiplying a suitable degeneracy factor to that of the 
scalar field. In the case of curved spacetime, however, 
the above statement is not valid in general. In this pa- 
per, we provide a complete and straightforward deriva- 
tion, using the Green function approach, of these degen- 
eracy factors, which are —4 and 15/4 for the untwisted 
and twisted BC, respectively, as expected. For a massive 
fermion, the Casimir energy density is plotted as a func- 
tion of the position y in Figj2j The asymmetry aspect 
of the surface divergence comes from the nature of the 
parity transformation in the fifth dimension, which can 
also be found in the flat case. 

To cure the surface divergence, an approach similar to 
that of 29J must be implemented where the surface terms 
are included. Under proper renormalization of the sur- 
face action term, the SUSY cancellation is not threatened 
by the power difference of the surface divergence between 
the scalar and the spinor fields. On the other hand, this 
infinity might be originated from the unphysical nature 



of the boundary condition. This may suggest the neces- 
sity of including the finite thickness (of string length) of 
the 3-brane in the treatment. In fact, the finite thickness 
effect of a de Sitter brane has been discussed in [3"3T - 
135] - Similar situation happens in the perfect conductor 
boundary condition for electromagnetic field [28[ [36] [37] , 
in which a strictly zero thickness boundary gives rise to 
a non-integrable divergence. For an imperfect conductor, 
such as a dielectric material, or a conductor with finite 
thickness, waves of sufficiently high frequency would pen- 
etrate into the material so that the precise location of the 
boundary would lose its meaning. Therefore, in reality, 
the divergence does not appear since the expression of the 
integrand is not universally applicable for all wave num- 
bers. The finite thickness of the 3-brane might provide a 
similar remedy to this divergence problem. 

Last, but not the least, there is the issue of self- 
consistency for the RS metric to retain its solution under 
the quantum corrected Einstein equations. As Knapman 
and Toms [5J commented in the case of the massless, 
conformally-coupled scalar field, different components of 
the stress energy-momentum tensor would in general in- 
duce different corrections to fc, if indeed an RS-type so- 
lution can be found. It was argued [B] that such correc- 
tion is exponentially small and therefore the RS solution 
remains approximately self-consistent. We found, how- 
ever, that in both cases of the massless Dirac spinor field 
and the conformally-coupled scalar field, the VEV of the 
stress-energy tensor takes the form (see (48) and (49 1): 



(T^) = 0(1) x e 3k\y\-5kKR k ^ 
(T yy ) = 0(1) x e 



5k\y \ — bkitR ^5 



(71) 
(72) 



where 0(1) is some quantities of the order unity. On the 
other hand, the cosmological constant term in the Ein- 
stein equation is AgMN ~ —M 3 k 2 gMN [2], where <?mjv 
is given by Q, so the ratio of stress-energy tensor to the 
classical one is about 



(k/M) 



3 p 5k(\y\--xR) 



(73) 



If k ~ M ~ Mpi , due to the requirement of no hierarchy 
problem, then the stress-energy tensor is comparable to 
the cosmological constant term in the region where y — ¥ 
irR, even though it is small with respect to k 5 in most 
regions. As a result, while the quantum correction of 
the metric is exponentially small at the UV-brane, that 
at the visible brane is of the same order as that of the 
classical Einstein equation with RS geometry. Thus we 
cannot treat the VEV of the stress-energy tensor as a 
perturbation to the semi-classical Einstein equation, and 
the metric of the form Q will not be a solution to the 
quantum corrected Einstein equation, not even in the 
approximate sense as argued in [B]. The self-consistency 
of the RS model that includes bulk fields is therefore in 
doubt. 

However [38], if k is slightly smaller than M as com- 
mented in [2], then the suppression (k/M) 3 would be 
rather significant even near the visible brane. Therefore, 



an approximate self-consistency of the RS-type metric 
may still be satisfied. 
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